CoUisional Effects in Isovector Response Function of Nuclear 

Matter at Finite Temperature 

S. Ayik\ A. GokaV, O. Yilmaz^, K. Bozkurt^ 
^Physics Department, Tennessee Technological University, Cookeville, TN 38505, USA 
'^Physics Department, Middle East Technical University, 06531 Ankara, Turkey 

(February 8, 2008) 



Abstract 

The dipole response function of nuclear matter at zero and finite temperatures 
is investigated by employing the linearized version of the extended TDHF 
theory with a non-Markovian binary collision term. Calculations are carried 
out for nuclear dipole vibrations by employing the Steinwedel-Jensen model 
and compared with experimental results for ^"^^Sn and ^^^Pb. 
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Typeset using REVT^ 



Giant resonances, in particular giant dipole resonances (GDR), in medium- weight and 
heavy nuclei have been the subjects of extensive experimental and theoretical studies during 
the last several years A large amount of experimental information is now available 
concerning the properties of GDR built on the ground and the excited states of the nuclei 
revealing the properties of the collective motion of nuclear many-body systems at zero and 
finite temperatures. The mean resonance energy is observed not to change much with the 
excitation energy, or the temperature, but the recent experimental investigations show that 
the width of the resonance becomes broader as excitation energy increases with a possible 
saturation at high temperatures. This temperature dependence of the GDR width is still 
one of the open problems in the studies of nuclear collective response and its damping 
mechanisms at zero and finite temperatures 0,^,01 • 

The theoretical investigations of the nuclear collective response employing the random 
phase approximation (RPA) theory have been quite successful in describing the mean reso- 
nance energies |p. However, the RPA theory is not suitable for describing the damping of 
collective excitations and therefore investigations based on the RPA theory have not been 
able to explain the increase of the width of GDR with temperature |^ . 

There are different mechanisms involved in the damping of the nuclear collective state. 
A part of the damping is due to the coupling of the collective mode to external degrees 
of freedom resulting in the cooling of the system by particle emission giving rise to escape 
width. Furthermore, the collective mode also acquires an intrinsic width as a consequence 
of its coupling to the internal degrees of freedom. The resulting spreading width which is 
thus due to the mixing of the collective mode with more complicated doorway states makes 
essentially the large part of the contribution to damping in medium-weight and heavy nu- 
clei. There are essentially three different theoretical approaches for the calculation of the 
spreading widths. In the first case, the temperature dependence of the width is explained 
by the coherent mechanism due to adiabatic coupling of the giant resonance with ther- 
mal surface deformations which is particularly important at low temperature. In the 
second approach, the mechanism of damping is due to the coupling with incoherent two 
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particle-two hole (2p-21i) states resulting from the increasing rate of collisions between the 
nucleons with temperature which is usually referred to as the collisional damping P,p!0 



The collisional damping is relatively weak at low temperature, but its magnitude becomes 
large with increasing temperature. The last mechanism is the Landau damping which is 
due to the spreading of the collective mode on non-collective particle-hole (p-h) excitations. 
Most investigations of nuclear response that have been carried out so far are based on the 
coherent damping or the collisional damping mechanisms [lTll , |T2 |. 



In this work, we perform a linear response treatment of the nuclear collective mode by 
including the collisional damping. The small amplitude limit of the extended the time- 
dependent Hartree-Fock (TDHF) theory, in which the collisional damping due to the in- 
coherent 2p-2h decay is included in the form of a non-Markovian collision term, provides 
an appropriate framework for investigating the damping widths of collective modes at zero 



and finite temperatures |T0| , pT| , p!2| , p!3| . We employ the extended TDHF theory to study the 
isovector response function with collisional effects of nuclear matter at zero and finite tem- 
peratures in semiclassical approximation using a simplified effective Skyrme force. 

The equation of motion of the single particle density matrix p{t) in the extended TDHF 
approximation is given by a transport equation 

th^^p-[hip),p]=Kip) , (1) 

where h{p) is an effective mean-field Hamiltonian and the right-hand-side represents a non- 
Markovian collision term, which can be expressed in terms of the correlated part of the 
two-particle density matrix as K{p) = Tr2[v,Ci2] with the effective residual interaction v. 
The correlated part of the two-particle density matrix Cu = P12—P1P2 where pip2 represents 
the antisymmetrized product of the single-particle density matrices, is given by the second 
equation of the BBGKY hierarchy. In the extended TDHF theory, the hierarchy is truncated 
at the second level by retaining only the lowest-order terms in the residual interactions, thus 
neglecting three-body correlations. Hence, the correlated part of the two-particle density 
matrix C12 satisfies the equation 



th^^C,2-[h{p),Ci2]=F,2 



(2) 



where the source term is given by 



Fi2{p) = (1 - - P2) V P1P2 - P1P2 V {1- pi)(l - P2) 



(3) 



In order to study the isovector collective response of the system, we include an external 
perturbation F{f,t) into the equation of motion, 



F{r,t)^TsF{r) 



(4) 



where T3 is the third component of the isospin operator and the frequency of the one- 
body harmonic perturbation operator contains a small imaginary part with the prescription 
CO ^ CO + irj in accordance with the adiabatic hypothesis. We obtain a description for 
small density fluctuations Sp{t) = p{t) — po in linear response treatment by linearizing the 
extended TDHF theory around a finite temperature equilibrium state density po, and this 
way we obtain 



ih^^Sp - [ho, 5p] - [5h + F(r, t),po] = Tr2[v, SC^] 



(5) 



where 5h — {dU /dp)Q5p represents small deviations in the effective mean-field potential. 
Moreover, the small deviation of the two-body correlations 5Ci2{t) satisfies 



i^^5Ci2 - [5h + F(f, t), C?2] - [/lo, <^Ci2] = 5Fi2 



(6) 



We look for a solution of Eq. 5 and Eq. 6 of the form 5p{t) — 5p{co)e~'^'^^ + h.c. where now 
the small density fiuctuations are given in terms of the proton and neutron density matrices 
as 5p{t) = Pp{t) — Pn(t)- We note that in the small amplitude limit, the collision term is also 
harmonic, 5K — tr2[v, SC12] — SK{co)e~'^'^^ + h.c, and in momentum representation we then 
obtain 



k k 
hu; + e(p- ■^)-e{p+-) 



k k 



fip- |) - f{p+ \) 



X <{ < p + ^ I 5/. I p - ^ > +2 < p+ ^ I F{r) I p - ^ > U< p+ ^ I 5K{uj) \ p~ ^ >, (7) 
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where f(q) — 1/ |l + e^[^*^^~^]| is the Fermi-Dirac occupation factor. For simphfied Skyrme 
interaction, the density fluctuation Sn{f, t) induces local changes in the mean fleld potential, 
therefore 



<p+-\6h\p- ->=2Vo<p+-\ 5p{w) \ p--> 



(8) 



where in terms of local proton and neutron mean field potentials 5h is expressed as 5h 
Up{r,t) — Un{r,t) — 2VQ5n{r,t). Moreover, since 



(9) 



and 



(10) 



we finally obtain 



Sn{k,w)- VoSn{k,w) + F{k) Ui{k,w) ^ VoSn{k,w) + F{k) U2{k,w 



(11) 



The function Ili{k,w), which is known as the unperturbed Lindhard function, is given by 



d^p 



f{p- I) - f{P+ I 



— * 

and the function n2(A;, w) is obtained from 

< p + I I 6K{w) I p - I > 



(12) 



[VoSn{k, w) + F{k)]U,{k, w) = I d'p 



hw - e{p+ I) + e{p- I) + ir) 



(13) 



The retarded response function which is defined by 



Sn{k,w) ^UR{k,w)F{k) 



(14) 



is then obtained as 



UR{k,w) 



l-VoIlo{k,w) 



(15) 



with Ilo{k, w) — li-iik, w) + Tl2{k, w) 
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From Eq. 6, it is possible to obtain a closed form expression for small deviation of 
two-body correlations (5Ci2, that is valid for collective vibrations (for details please refer to 



i rt 

n 



dt'plple~''''''^'-''\5^{t'),vy''''^'-''\l - p?)(l - pI) + h.c. 



(16) 



where 5^{t) is the distortion function associated with the single-particle density matrix, 
and it is related to the small vibrations in the single-particle density matrix 5p{t) according 
to 5p{t) = [5$(t),po]- We then obtain the expression for the linearized collision term by 
evaluating the matrix element < p + | | SK{u) | p — | > in which we retain fc-dependence 
only in distortion function Then, the coUisional response function 112 (/c,w) can be 

expressed as 



-3 ,3 -3 ,3 f^QV W{12; 34) /1/2/3/4 - /1/2/3/4 

/o *A3 / ^ p^^ ^ i ^ — ^ — r~ ^'^ 

[2nhy J V 2 / vr kw — — + ei + €2 + ir] 



where = 1 - fi, AQ = Qi + Q2 - Q3 - Q4 with Qi = 1/ \hw - e{pi + |) + e{pi - |) 
and W(12;34) denotes the basic two-body transition rate, which can be expressed in terms 
of the spin averaged proton-neutron scattering cross section as 

1 4h 



Vr(12;34) 



(27r^)3 m2 \dnj 



^^1 5^(pi +P2 -P3 -P4) • 



(18) 



pn 



The strength distribution function is obtained from the imaginary part of the retarded 



response function |I4 



S{k, w) 



TT 



-ImIlR{k,w) . 



In our calculations, we employ a simplified Skyrme interaction 

V = toil + xoPa)6{r) + ^^3(1 + X3P.)p"(i?)5(f) 



with f = f{ — 7^2 and R = {f\+ 'i%)/2. The local potential for protons is then given by 
f/p(f, t) = to (^1 + ^xo^ p(f, t) - ^0 Q + pp{r, t) 

(2 + a) {l + ]^x^pp{r,t) 



(19) 



(20) 



+^i3P"(r,t) 



-2 ( ^ + p{f, t) - a Q + X3 



p(f, t) 



(21) 



with a similar expression for neutrons. In linear response approximation, the coupling 
constant Vq for dipole vibrations becomes 

where po is the saturation density of nuclear matter. In our analysis we consider in 
particular the Skyrme SLy4 force with the parameters |15| to = —2488.91 MeV.fm^, 



ts = 13777 MeV.fm'^''^, xq = 0.834, = 1.354 and a = 1/6, which results for Vq in 
the value Vq = 85 MeV. 

In order to apply our results to finite nuclei, we work within the framework of Steinwedel 
and Jensen model for nuclear dipole oscillations 0. In this model neutrons and protons 
oscillate inside a sphere of radius R given by the expression 

Pp{f,t)-pn{f,t)=Fsm{k-r)e'^' , (23) 

the total density remaining equal to the saturation density po of nuclear matter and the 
wavenumber k is given by A; = tt/2R. We apply Steinwedel and Jensen model to GDR 



m 



120 



Sn and "^^^Pb, and we take R = 5.6 fmk = 0.28 fm ^ for ^"^^Sn and R = 6.7 fm 



k = 0.23 /m-i for ^°^Pb according to R = 1.13 A^^^. 

As a result of the approximate treatment, the collisional response function Il2{k,uj) 
has a singular behavior arising from the pole of the distortion functions, Qi = 
1/ hw — e{pi + |) + e(pi — |) . We avoid this singular behavior by incorporating a pole 
approximation. In the distortion functions, we make the replacement uj ^ ud — i'^/^ where 
uj£) and r are determined from 1 — Voni(A;,co') = at each temperature that is consid- 
ered. Furthermore, we neglect the real part Rell2 of the function 112 in our calculations. In 
the calculations of the collisional response function, using conservation laws and symmetry 
properties, it is possible to reduce the twelve dimensional integrals to five fold integrals by 
incorporating the transformations into the total momenta P = pi + P2, P' = P3 + Pa, and 
relative momenta q = {pi —p2)/2, q' = {p^ —p^/2 before and after the collisions. The inte- 
gral over P' can be performed immediately. The delta function 5{huj — e' + e) in Iml\.2{k, w) 



where e = q'^/m and e' = q' ^/m are the energies of two particle system in the center of 
mass frame before and after the colhsion makes it possible the reduce the integrals further 
using familiar methods from the Fermi liquid theory ||T^. Then, we evaluate the remaining 
five dimensional integrals numerically by employing a fast algorithm. In the evaluation of 
momentum integrals, we neglect the angular anisotropy of the cross sections and make the 
replacement {da/dVt)^^ (Tp„/47r with (Tp„ = 40 mb. 

We show our results for the response function with and without the collision term in 
Fig. 1 for ^'^^Sn and in Fig. 2 for "^^^Ph as a function of experimental temperature T* 
where we also present the comparison with the normalized experimental data 0]. The 
experimental temperature T* is related to the temperature parameter in the Fermi-Dirac 
function /(e, T) as T = T* ^ ap, where ap denotes the energy dependent empirical level 
density parameter and ap = ATc^/Atp denotes Fermi gas level density parameter P,p!6|. 
In our calculations, we use the temperature values T in the Fermi-Dirac function that are 
related to the experimental temperatures in this manner. 

We first note that, at the RPA level that is without the contribution of the collision term 
in calculation of the nuclear response, the position of the peak of the response functions do 
not change with temperature. As a matter of fact, for ^osp^ ^j^g peak is at = 12.3 MeV 
for T* = 1.34, 1.62, 1.85, 2.05 MeV while for ^^o^^ jt occurs at u = 14.3 MeV for 
T* = 1.36, 2.13 MeV and at u = 14.2 MeV for T* = 2.67, 3.1 MeV. This behavior of 
the peak energy with temperature is in accordance with the experimental results where it 
is observed that the mean-energy of the dipole response is almost constant for ^osp^ when 
T* changes between 1.3 to 2.0 MeV while a decrease of 1.5 MeV is observed in ^"^^Sn when 
T* changes from 1.2 to 3.1 MeV ||^. However, the average positions of the peak values of 
the strength functions are slightly below the experimental values, which are u = 15.4 MeV 
for ^'^^Sn and uj = 13.4 MeV for ^^^Pb. This discrepancy may be due to the nature of the 
effective Skyrme force that we employ. Moreover, the value for k = tt/2R that is used in 
Steinwedel and Jensen model depends on the value of Rq used in _R = RqA^/"^., and changing 
k somewhat also produces a change in the position of the peak but in general the above 
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conclusions are not affected. Furthermore, These results for the peak position of the strength 
functions are in accordance with the earlier RPA calculations. 

Dotted lines in figures 1-2 show the strength functions including coUisional damping 
mechanism. Since, we neglect the real part of the coUisional response the peak values of the 
strength do not change, but the coUisional mechanism introduces a spread, in particular at 
the high frequency side of the strength functions and this spread becomes more pronounced 
with increasing temperature. In order to illustrate the effect of the collision term more 
drastically, in figure 3, we show the response function with and without the collision term 
for ^"^^Sn and for 208p^ temperatures T=2,4 MeV. 

Rather simple description presented in this paper is able to explain certain aspects of 
giant dipole excitations in ^^^Sn and "^^^Pb, but do not produce a good description of the 
experimental strength functions as a function of temperature. One important element miss- 
ing in the calculations is the coherent damping mechanism due to coupling dipole vibrations 
with low frequency collective surface modes |l8| , p!9| , p0| . This mechanism is especially im- 
portant for describing the details of the strength distributions at low temperature. In the 
continuation of this work, we plan to improve this simple description by incorporating the 
coherent damping mechanism in Thomas-Fermi approximation. 
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FIGURES 




FIG. 1. The GDR strength function of Sn. Sohd and dashed hnes show the response function 
without and with the colHsion term, respectively. The normahzed data is taken from [4] . 
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FIG. 2. The GDR strength function of °°Pb. Sohd and dashed lines show the response function 
without and with the coUision term, respectively. The normalized data is taken from [4] . 
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FIG. 3. The GDR strength function of ^^OSn and ^ospb at T* = 2, 4 MeV. Sohd and dashed 
hnes show the response function without and with the colhsion term, respectively. 
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